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JUNCTION CAPACITANCE .AS A FUNCTION OF VOLTAGE FOR DIFFUSED 


p-n DIODES WITH EXPONENTIAL DOPING GRADIENTS 
by Ira T. Myers and Julian F. Been 
Lewis Research Center 

SUMMARY 

A method was developed for calculating the capacitance of a diffused p-n junction as 
a function of voltage that can be used in place of the more complicated methods now 
available. A single universal curve, based on an exponential doping gradient, permits 
the calculation of the capacitance of junctions with a wide variety of base doping levels 
and diffusion doping gradients. In particular, the method permits the calculation of ca- 
pacitances at intermediate voltages where the junction behaves neither as a linear graded 
junction nor as an abrupt junction, but is intermediate between the two cases. Calculated 
and measured capacitances as a function of voltage were compared for a silicon power 
diode with good agreement. 


INTRODUCTION 

The capacitances of linearly graded and abrupt junctions are treated in many text- 
books (refs. 1 and 2). A reverse biased diffused p-n junction will normally show a 
(V - dependence of capacitance at low voltages and a dependence at high 

voltages . The "built-in” or diffusion voltage of the junction is Vg. 

Between the low and high voltages a transition region exists where the capacitance 
fits neither the linear graded (V - nor the abrupt case. Lawrence and 

Warner (ref. 3) calculated this case using the complementary error function and the 
gaussian distribution as doping distributions. This report describes a similar calculation 
using an exponential doping distribution. Justification is given for using the exponential 
doping distribution results for all diodes, no matter what the doping distribution actually 
is. The advantages of using the exponential doping results, as given in this report, are 
shown to be improved clarity of interpretation of the doping gradient, and a considerably 
more straightforward calculation procedure. This report concludes with a comparison 



of calculated and measured capacitances as a function of voltage for a silicon power 
diode (ref. 4). 


CALCULATION OF CAPACITANCE AGAINST VOLTAGE 
FOR A DIFFUSED p-n JUNCTION 

Consider a junction where the dopant has been added such that the added impurity 
follows an exponential distribution. For example, take p type base material, with a 
concentration of N^ dopant atoms per cubic centimeter. The concentration N^ of the 
n type dopant atoms is given by 


“d = V”’' (« 

where N„ is the surface concentration, X the distance from the surface, and b a con- 

s 

stant. The net doping Nj^ - N^ is given by 

N = Nge-*^^ - Nq (2) 


At the junction N = 0 and X = X^ so that 


No 


-bX, 




(see fig. 1). The measured capacitance C of the junction is 

C = ^ 
dV 


(3) 


(4) 


where V is measured in the reverse direction. 


The capacitance 


C is also given by 


C = 


eCoA 


W 


(5) 


where 


W = Xg - Xj 


( 6 ) 


2 



and where W is the depletion region width, ee^ the permittivity, and A the area of the 
junction. 

With applied voltage the depletion region expands an amount dXj in the n region 
(X^ -*■ Xj - dXj^) and an amount dXg in the p region (X 2 — X 2 + dX 2 ) . The amount of 
positive charge uncovered by the expansion of the depletion region in the n region is 

dQ^ = qA (dX^) (7) 

Similarly, the amount of charge uncovered by the expansion of the depletion region in the 
p type base is 



= qA 




dX2 


In order to maintain charge neutrality over the entire depletion region, 


dQ„ = dQp 


( 8 ) 


(9) 


The total dW is given by 


dW = |dX2l + |dX^| (10) 

Next, equations (4), (5), and (7) to (10) are used to eliminate dQ and dQ and to ex- 
press dW in terms of dV: 


dW 



( 11 ) 


This is a key equation in the development. It will be integrated to give W and therefore 
C as a function of V. 

Before equation (11) can be integrated, Xj and X 2 must be expressed in terms of 
W. The quantity X 2 is easily expressed in terms of X^ and W as 

Xg = W + X^ (12) 
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so that it is only necessary to find the relation between X| and W. This is obtained by 
equating the positive charge on the n side of the junction to the negative charge on the 
p side of the junction: 


qA 



N dX = qA 




(13) 


After integration and simplification the expression for is 

-bX. bN W 
e 2 

(1 - e-'^X 


(14) 


Substituting equations ( 14 ) and (12) into equation ( 11 ) gives 


dW = 




dV 


qN^W 


1 - e 


-bW 


1 - e 


-bW 


1 - e"^^ - bWe"^’^ bW 


+ e 


-bW 


(15) 


where dV is positive for increasing reverse voltage. 

As stated previously, W and V are variables of integration. Introducing dummy 
variables of integration (fi for bW, and for ¥) allows the use of W and V for in- 
tegration limits. After some simplification, 


^2(1 + e'^) 
1 - e'^ 


(1 - e"^) 


- 1 dJ2 = 


b^ee^ 

In" 


d(p 


( 16 ) 


This expression was integrated over the depletion region. The limits for f2 were 
= 0 and S2 = bW, where W is the depletion region width at a voltage V. The cor- 
responding limits for <p are - Vg (corresponding to W = 0) and the applied voltage V. 
Therefore, 



n(l + e~^) _ f2^e'^ 

(1 - e-") (j . 


dS2 = 


b^e€j. 



(17) 
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where the substitution z = bW has been made. This may be easily Integrated by inte- 
grating the second term in the first integral, 



(1 - 


dn 


by parts and by combining the results with the other terms prior to integrating them. 
The result is 


1 - e-^ 2 qN„ 


( 18 ) 


Define the function G(z) as 


G(z) = 



1 - e 


-z 



( 19 ) 


where z = bw as defined previously. The function G(z) is always positive. Thus 


G(z) = 


b^ee^CV + Vg) 


qN, 


( 20 ) 


Then 


z = bw = 


G 


■lA ^S(v-^Vb )^ 

qN. 


( 21 ) 


_i 

where G is the inverse function of G. It was not possible to find the inverse of G 
and express it in closed form algebraically (i. e. , an expression f(G) = z, which would be 
the solution for z of ec^. (19)). Instead , G was plotted against z (fig . 2) and tabulated 
(see table I) so that for any given value of G the value of z may be determined. 

In the range 0. 1 < G < 10, the quantity z (or bw) can be represented, within ±5 per- 
cent, as 
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z = 1.80 + 0.60 


( 22 ) 


where G = b^eeQ(V + Vg/qN^. For G < 0. 1, 

z = (12G)^/^ (23) 

is a close approximation, while for G > 10, 

z = 1 + Vl + 2G (24) 


can be used. 

The capacitance against voltage can now be directly calculated: 

€€ Ab 
C = -^ — 

z 


(25) 


where the quantity z contains the voltage dependence. For the limiting case of low 
voltages, where z = (12G)^^^, the capacitance reduces to the expected cube root de- 
pendence 


C 

A 


12(V -f Vg) 


1/3 


Similarly at high reverse voltages, where z = 1 -i- y/l + 2G \/^, 


C 

A 


ee^qN 


0 


1/2 


2(V -f Vg) 


(26) 


(27) 


which is the expected square root dependence. 

In summary, to calculate the capacitance at a given voltage V, first calculate 
b\eQ(V + Vg)/qNQ, then find z from figure 2 or table I corresponding to it, and finally 
calculate C from 


C = 


ee^bA 


z 


(28) 
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If Vg is not known, it can be calculated from the work of Nuyts and Van Overstraeten 
(ref. 4). A fit to their results gives 

Vg = 0. 080 logj^QbN^ - 1. 08 (29) 

17 23 4 

in the range of bN^ between 10 and 10 (bN^ is the doping gradient in cm" at the 

junction) . 


SAMPLE CALCULATION 

To illustrate the method, including the units to be used, a sample calculation is 
given. Consider a diode with the following characteristics: 

2 

A junction area, 0. 1 cm 

V applied voltage, 3.0 V 

4 -1 

b logarithmic doping gradient, 10 cm 

15 3 

Nq base doping level, 10'‘’Vcm 

19 -4 

UNq doping gradient at junction, 10 cm 

Since Vg is not known, calculate it from the Nuyts and Van Overstraeten expression 

Vg = 0. 08 logjQ bN^ - 1. 08 = 0. 44 V 

Then 

G = 06x10~^^)(3. 0 + 0. 44) ^ 2 28 

1.602x10"^^ (10^^) 

and z = 3. 2 from figure 2. From this value of z 

c = zz 330x10"^^ F or 330 pF 

z 3. 2 

COMPARISON OF CALCULATED AND MEASURED RESULTS 
The technique just described was used to calculate the capacitance as a function of 
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voltage for a p-i-n type diode number SIN1189 (ref. 5). This was a 35-ampere diode 

with a junction area of 0. 455 square centimeter and a base doping level of about 1. 2x10^^ 

3 

cubic centimeters. The value of b was about 5x10 per centimeter. 

For the comparison, was determined from experimental measurements of ca- 
pacitance by fitting at the high voltage end of the curve, and b was determined by fitting 

1 4 

at the low voltage end. The N„ as determined for the diode used was 1. 1x10 per 

o 3 

cubic centimeter and b was 5. 1x10 per centimeter. The diffusion voltage V„ was 

3 ^ 

determined both by extrapolation of a 1/C plot against voltage, and from b using 
Nuyts and Van Overstraeten's work. The value of Vr, as determined graphically from 
the l/C plot was 0. 35 volt; the calculation from Nuyts and Van Overstraeten's work 
gave 0. 34 volt. 

Figure 3 gives a comparison between calculated values and measured values of ca- 
pacitance. The agreement is good. 

Table n gives a comparison of experimental results with the capacitances as calcu- 
lated using the three different types of doping gradients (exponential, from this report, 
and complimentary error function and gaussian from Lawrence and Warner’s work 
(ref. 3)). To compare the goodness of fit for the three doping gradients used, the ratio 
^calculated/ ^measured averaged over the voltage points from 1 to 100 volts. For 
a perfect fit this ratio should be 1. 00. The value of this ratio, and its error, give an 
estimate of how the model agrees with experiment (see table HI). Ninety percent con- 
fidence limits are used on the errors. 

The curves were fitted at zero bias and at 200 volts in each case. There is less 
difference between the calculations using different types of gradients than between cal- 
culated and experimental values. The particular kind of doping gradient assumed does 
not make much difference in the calculated capacitance. The variation in calculated ca- 
pacitance from one doping gradient type to another is within the calculated errors. This 
could be expected, since at very low biases, with a narrow depletion region, any gradient 
can be fitted with a straight line. On the other hand, at high voltages, nearly the entire 
depletion region is in the base, so that again the type of gradient is relatively unimpor- 
tant. Therefore, the simplest to use gradient, namely the exponential, may be used 
under a wide range of conditions . 


SUMMARY AND CONCLUSIONS 

A method of calculating capacitance as a function of voltage for an exponentially 
doped p-n junction is described. The method is significantly easier to use than pre- 
sently available methods where a complementary error function or gaussian doping gra- 
dient is assumed. Comparison with a power diode gave good agreement with experi- 
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imental results. Calculated capacitances using exponential, complementary error func- 
tion, and gaussian distribution all gave results in agreement with measured values 
within the calculated errors. This indicates that the dopant distribution assumed is not 
critical, and that the more convenient and easier to interpret exponential distribution can 
be used for capacitance against voltage calculations . 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 7, 1970, 

120-60. 
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TABLE I. - NUMERICAL VALUES OF G{z) FOR VALUES OF z FROM 0 TO 9.9 


z; 

0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 j 

0 

a 

8. 3319X10'^ 

6.6622x10'^ 

2.2466x10'^ 

5.3192X10“® 

1.037X10’^ 

1. 7893x10*® 

2.8353x10*® 

4.2218x10*® 

5.9945x10*® 

1 

8. 1977x10'^ 

1.0874x10"^ 

1. 4066X10*' 

1.7813x10*' 

2.2153X10*' 

2.714X10*' 

3.2760x10“' 

3.9096x10*' 

4.6163x10*' 

5. 3990X10"' 

2 

6.2607x10'^ 

7. 2039x10’' 

8.2311X10*' 

9.3447X10*' 

1.054T 

1.1839 

1.3224 

1.4702 

1.6276 

1. 7947 

3 

1*9716 

2.1583 

2.3551 

2.5620 

2. 7791 

3. 0064 

3.2441 

3.4920 

3. 7504 

4.0192 1 

4 

4.2985 

4.5883 

4. 8885 

5. 1993 

5. 5206 

5.8525 

6.1949 

6.5478 

6.9112 

7.2851 

5 

7. 6696 

8.0645 

8.4670 

8. 8859 

9.3123 

9. 7491 

10. 196 

10. 654 

11. 122 

11.601 

6 

12.089 

12. 589 

13.098 

13.618 

14. 148 

14.689 

15.239 

15. 800 

16. 372 

16.953 

7 

17. 545 

18. 147 

18. 759 

19.381 

20.013 


21.309 

21. 972 

22. 645 

2S. 328 ^ 

8 

24.021 

24. 725 

25.438 

26. 162 

26.896 

27.640 ! 

28. 394 ! 

29.158 

29.932 

30. 716 

9 

1 31. 510 

32.314 

33.129 ! 

33.953 

34.787 

35.632 

36.486 i 

37, 351 

38.225 1 

39. 110 













table n. - COMPAKBION OF THREE 


METHODS OF CALCULATION OF CA- 


PACITANCE WITH EXPERIMENT 


Applied 

Calculated capacitances 

Experimental 

voltage 

^erfc 

C 

gauss 

c 

expon 

capacitance , 
^exper 

0 

1300 

1300 

1300 

1300 

1 

840 

770 

790 

840 

2 

640 

640 

645 

680 

5 

480 

470 

466 

490 

10 

380 

360 

357 

375 

20 

270 

270 

270 

290 

50 

190 

180 

181 

190 

100 

137 

136 

131 

132 

200 

95 

95 

95 

95 


TABLE m. - RATIO OF CALCULATED TO MEASURED 


CAPACITANCES FOR DIFFERENT DOPING GRADIENTS 
























Figure 2» - Inverse junction capacitance function 


z versus the applied voltage function 


G. 
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